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Abstract 

We present a semiclassical nonlinear field equation for the con- 
fining field in 2+1-dimensional U{1) lattice gauge theory (compact 
QED). The equation is derived directly from the underlying micro- 
scopic quantum Hamiltonian by means of truncation. Its nonlineari- 
ties express the dynamic creation of magnetic monopole currents lead- 
ing to the confinement of the electric field between two static electric 
charges. We solve the equation numerically and show that it can be 
interpreted as a London relation in a dual superconductor. (UFTP 
preprint 405/1996) 



1 Introduction 

In the last few years, several studies of confinement in lattice gauge theories 
have succeded in observing not only the confinement potential |^, but 
also the formation of flux tubes between two static charges 0, ^ . The most 
promising mechanism for this effect is the dual superconductor hypothesis 
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1^, which assumes that the electric field is confined by magnetic monopole 
currents forming the walls of the fiux tube, just like superconducting elec- 
tric currents confine the magnetic field in an ordinary superconductor. The 
monopoles are presumed to originate dynamically from the periodicity of the 
Hamiltonian as tunneling effects between neighboring minima [|, Q . Abelian 
monopoles occur naturally in U{1) lattice gauge theory (compact QED). 
As U{1) is a subgroup of SU{2) and SU{'i), so-called "Abelian" monopoles 
(monopoles in the maximally Abelian projection) |jlO|] may also play a role 
in the QCD mechanism of confinement. The dynamic creation of magnetic 
monopoles and their description by a London relation has been observed in 
several studies [O, [l^, O, |l4l based on DeGrand and Toussaint's identi- 



fication prescription [O for monopoles. 

Recently, Zach et. al. demonstrated that confinement of the electric 



field inU{l) lattice gauge theory can be effectively described by the Maxwell- 
London equations of a dual superconductor. This model has only a single 
parameter, the London penetration length. It is, however, phenomenological 
in the sense that it is not derived from the underlying microscopical lattice 
system. 

Most studies have employed Wilson's Euclidean picture of lattice gauge 
theory [jl6|, simulating it as a statistical system using Monte Carlo methods. 
While this offers easy access to the confining potential (string tension) by 
measuring the expectation value of a Wilson loop, and to the monopole 
density by means of the DeGrand- Toussaint prescription, measurements of 
the fields involves a more complex correlation function (basically between a 
loop and a plaquette) and thus more statistical fiuctuations. 

An alternative approach is to make use of the Kogut-Susskind Hamilto- 
nian formulation [Q, considering the theory as a quantum system with many 
degrees of freedom and applying methods from many-body theory as devel- 



oped in e.g. computational chemistry or nuclear physics [19, EQ, Ell 23, 24 



In particular, the Hamiltonian formulation provides more easily for a semi- 
classical static solution than the Euclidean version. 

In the following, we perform an ah initio derivation of a semiclassical field 
equation in 2 + 1 dimensions that incorporates the prime feature of the gauge 
group U{1) — the periodicity of the action that leads to the dynamic creation 
of effective monopole currents. The resulting nonlinear field equation exhibits 
confinement of the electric field caused by magnetic monopole currents. It 
is microscopic in the sense that it is formulated in terms of the plaquette 
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variables (i.e. the vortex potential of the electric field) and does not make 
use of phenomenologically introduced monopole fields. 

We first review the formulation of classical electrodynamics on the lattice. 
After introducing the Kogut-Susskind Hamiltonian, we move to the plaquette 
representation (which only exists without constraints in 2+1 dimensions) 
and then perform a truncation of the Hilbert space to obtain a governing 
equation for plaquette wavefunctions. By studying single-plaquette solutions, 
we are led to an effective field equation in terms of the plaquette variables. 
Its linearized form is the London relation. The numerical solution shows 
confinement of the electric field. We discuss the consequences of this for the 
confinement of charges in lattice gauge theories. 



2 Derivation of the effective field equations 

2.1 Classical electrodynamics on the lattice 

We are interested in the distribution of the electric field Ei on the links / of 
a 2-dimensional rectangular lattice in the presence of static charges. In the 
continuum, the electric field is determined by two conditions 

V-E{r)^Pe{f) , (1) 
V X E(r} = . (2) 

The first is Gauss' law, expressing that the electric charge density is the 
source of the electric field, the second Stokes' law that there are no closed 
lines of force (when we consider the static case without currents). These 
two equations determine the electric field completely. If the field exhibits 
confinement, one or both of them have to be modified. As the first equation 
is linked to gauge invariance of the theory (which is not violated), the only 
possible modification is to replace the zero on the right-hand side of the curl 
equation by a quantity — J, interpreted as the density of magnetic monopole 
currents. 

On the lattice, the electric field is represented by its fiux Ei along the 
directed links I. Gauss' law then becomes a condition on the link fields Ei 
joined at a site s 

J2^s,iEi^pi (3) 
I 
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where es,i is +1 if the hnk / leads away from site s, —1 if it leads to site s, 
and zero otherwise. In other words, the difference of inbound and outbound 
flux at any site equals the charge at this site. 

Stokes' law turns into a condition on each plaquette p, 

J2e,^,Ei = -J, (4) 

I 

where ep^i is +1 or —1 if the directed link / counts positively or negatively, 
resp., when going counter-clockwise around the plaquette. This condition 
states that the electric flux around a plaquette is generated by magnetic 
monopole currents. 

2.2 Quantum electrodynamics on the lattice 

In quantum field theory, the electromagnetic field is governed by the Kogut- 
Susskind Hamiltonian for the gauge group U{1): 

« = ~ E ^ + A E (l-Ree«^) , (5) 

where A = l/g"^, 0/ is the group parameter associated with the link /, and 
Gp the counter-clockwise sum of the link parameters around the plaquette p, 

I 

The Hamiltonian operates on a Hilbert space spanned by basis states 

|0)= n -vr<ez<7r . (7) 

The first (kinetic) term in eq. (^) is associated with the electric field, the 
second (potential) term with the magnetic field. Without the nonlinearity 
introduced by the higher-order terms of the cosine in the magnetic field term, 
this Hamiltonian would describe the motion of (nonrelativistic) particles in 
a harmonic potential (noncompact QED). In compact QED, the harmonic 
potential is replaced by a periodic potential that can be effectively described, 
as we shall see, by the dynamic creation of monopole currents. 
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As we are interested in the electric field, it is advantageous to switch to 
eigenstates of the momentum operator 

given by 

(0,|n,) =e^"'®' . (9) 

The Hilbert space will then be spanned by basis states \n) whose label n 
stands for a set of integer numbers rii, one for each link I: 

\n)= n 1^')' ^^ez . (10) 

In this representation, the Hamiltonian can be written 

H=l E + A E i^-l^p-l^^) ' (11) 

/Glinks pgplaquettes 

where n represents the momentum operator (§) with the matrix element 

{n'\h\n) = n6n,n' , (12) 
and the group element operator Up 

Up = e'®^ (13) 

with the matrix elements 

{n'\Up\n) = 5n',n+p , {n'\U+\n) = 6n',n-p , (14) 

where n' — p represents the configuration n' with all links in the plaquette p 
raised or lowered by one unit according to whether they count positively or 
negatively when traversing the links of the plaquette counter-clockwise. 

Note that the periodicity of the action has now transformed into the 
discretization of the dynamic variable n. This is the only nonlinearity in this 
equation. 

The Hamiltonian commutes with the operators 

= E (^s,ini, s e sites (15) 

/G links 
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that generate gauge transformations at the lattice sites s. In the electric 
field representation, these operators are diagonal and their eigenvalues give 
the total electric field flux out of the sites s. By means of Gauss' law (|^), 
they are associated with static charges sitting at the lattice sites s. The 
eigenstates of the Hamiltonian separate into different sectors, each associated 
with a certain distribution of integer charges over the lattice sites. We are 
especially interested in the string sector with two opposite charges a certain 
distance away from each other. 

In the strong-coupling limit A — > 0, the Hamiltonian is dominated by the 
electric term. The lowest eigenstate in the string sector is the shortest path 
of raised links connecting the two charges. As the magnetic contribution is 
switched on by raising A, perturbations are created from the original ground 
state by the action of the plaquette raising and lowering operators Up and 
Up. 

2.3 Plaquette representation 

In two dimensions, self-duality of the lattice allows us to introduce new vari- 
ables rrip associated with plaquettes p to replace the link variables rii by the 
relation 

ni = Yj(^p,i^p + (16) 

p 

As there are exactly two plaquettes p for which e^^/ is nonzero, namely the 
plaquettes to the left {ep^i = 1) and to the right {ep^i = —1) of the link /, 
(|T6|) expresses each link as the difference of the two adjacent plaquettes. The 
corresponding continuum relation is 

E = V X Th + N (17) 

where the vector field rh is in the z-direction such that E is in the x-y- 

plane. Thus m is the vortex potential of the electric field. Gauss' law is 
automatically satisfied 

V ■ ^ = V ■ iV = (18) 

given an appropriate choice of A^ as a representative for the sector of the 
theory, e.g. a Dirac string for the two-charge sector. There is a residual 
gauge invariance under global shifts of the rrip. 
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The original Hilbert space is represented by the direct product of the set 
of all plaquette variable configurations rrip and one configuration A*"; represen- 
tative for each sector of the theory. It thus allows us to focus our attention 
on one sector by fixing A*"; while admitting all configurations rrip. 

Turning first to the classical equations @ and we note that, while 
Gauss' law is automatically satisfied, Stokes' law becomes 

- ^ rup' + Np = -AiTUp + Np = -Jp (19) 
(pp') 

where the sum runs over all neighbors p' of the plaquette p, Al is the lattice 
Laplace operator on an unit lattice, and Np is defined as the lattice curl of 

Np = T.^P,iNi ■ (20) 
I 

In particular, if Ni is a string connecting two static charges, A'p is +1 on the 
plaquettes just above the string and —1 just below the string. 

The quantum Hamiltonian in the plaquette representation reads 



H = 2j2mi-Y: ^mp, + ^ E + E ^p^p + ^ E (i " ^Up - 

V (pp') I V p \ ^ 

(21) 

where (pp') sums over neighboring plaquettes p' . We note that, while the 



diagonal (potential) part of (|2T|) is quadratic and thus amenable to a Gaus- 
sian' solution, the operators Up are highly nonlinear as they implicitly contain 
higher powers of the derivative operator. This, of course, is an expression of 
the fact that m is confined to integer values which, in turn, results from the 
periodicity of the original f/(l) action and will ultimately give rise to a mod- 
ification of Stokes' equation. We also note that this restriction is enforced by 
the definition of Up alone so that we may consider m a continuous variable. 

2.4 Truncation 

To solve the eigenvalue problem of eq. (0), we truncate the eigenvalue 
equation to states that are direct products of single-plaquette states, 

i^)=ni^p) • (22) 
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To find the governing equation for plaquette g, we project onto tlie state 

■ (23) 

Tfiis leads to tfie relation 

+ 



where the Hamiltonian was decomposed according to 

H^J2Hp+ E Hpp, + Ho (25) 
p pp',p^p' 

{Hq accounts for the constant part in the Hamiltonian). In this approxi- 
mation, we find as the governing equation for the plaquette wavefunction 

\i/jp) = Ep\i/;p) , 

(26) 

that is, a single-plaquette equation with the next-neighbor interaction re- 
placed by a mean-field term. The single-plaquette energy Ep is related to 
the total energy by 

Ep^E-Y^{H,)-^^{H,,,)-Ho . (27) 

q^p q^p q'^p 

Summing over all Ep would double-count the next-neighbor interactions. The 
correct ground-state energy is given by the expectation value of the Hamil- 
tonian. However, since we have dropped all correlations between plaquettes, 
its value in this approximation cannot be trusted. 

This approach can be extended to basis functions defined on a blocked 
lattice, as would be required for a renormalization group analysis. 
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E i^plHpl'ipp) + E E {'^p'^p'\Hpp'\lpp1pp') I \lpq/ 

\p¥=g p^<ip'¥=q 

Hql'^q) + E {^p\HqP + Hpq\lppiljq) 



p¥=g 



E\^q 



(24) 




2.5 Single-plaquette solutions 

We rewrite the single-plaquette equation 



H,,p, = 27h' - Jm + Xl^l- -U - -U^j (28) 

introducing a generic source term —Jm that accounts for both the original 
source term and for the next-neighbor interaction: 

Jp= i'^p') -^P = ^LTUp + Arup - Np . (29) 
(pp') 

We shall discuss this equation by considering its solution in the strong- and 
the weak-coupling limit. 

In the strong-coupling limit A ^ 0, the Hamiltonian reduces to 

2m^ - Jm . (30) 

As it is diagonal in m, eigenfunctions are concentrated on a single value of 
m: 

{■m\ljjn) = Sn,m ■ (31) 

The state associated to the lowest eigenvalue is found by minimizing the 
energy over the integer numbers: 



-Emin = min„ez(2m - Jm) ^ — — + 



J 1 

4 + 2. 



4^ (32) 



where [x] denotes the smallest integer not larger than x and thus Jo the 
multiple of 4 closest to J. 

In the weak-coupling limit A ^ oo, we assume that the ground-state wave 
function {m\ip) is spread out sufficiently wide that it can be approximated 
by a Gaussian: 

(m|z^) ^ iVexp |^_ (^ (33) 

with a normalization constant N, first moment mo and second moment a. 
We have 

{m\U\^) = {m + = e-(^+2("*-"*°»/^'^(m|7/;) (34) 
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and thus 



{m\H\ 



+ A 



2rn? — Jm 
1 
2 

E{m\il) 



\ (,-^1^'^ |^g(m-mo)/2o- _|_ g-(m-mo)/2o- 



(35) 



Assuming that a ^ 1 and considering values of m in the vicinity of mo, we 
can assume \{m — mo)/cr| ^ 1 and expand the exponentials. Keeping terms 
up to l/o"^, we find the eigenvalue equation 



2m — Jm + A 



8(T 



' (m - mp)^ 
8^2 



To satisfy this equation for any m, we set 



TA 



a 



mo 



J 



(36) 



(37) 



The eigenvalue is 



J2 



(38) 



We thus find both in the strong- and weak-coupling limit that the expec- 
tation value of the plaquette variable m is approximately 



m 



= (m) ^ — 



(39) 



However, while this becomes exact in the weak-coupling limit, in strong cou- 
pling the wave function feels the restriction of m to integers. Fig. |l] shows 
the deviation of the numerically calculated expectation value of m from the 
weak-coupling expression as J is varied. It exhibits a sawtooth shape in 
strong-coupling that is softened to a sine-like shape in the intermediate range. 
The cause of the deviation is the restriction of m to integer values which is 
a direct consequence of the periodicity of the action. 

We can thus parametrize the nonlinearity of the single-plaquette Hamil- 
tonian by writing 

J . /TT 



(mapprox.) = - - asin l^-J 



(40) 
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J 



Figure 1: Deviation of the expectation value of the plaquette variable (m) 
from its weak-coupling limit as a function of the external source J at different 
couplings A. This quantity is a measure for the nonlinearity of the system. 



The parameter a is a function of A and vanishes in the weak-coupling limit. 
In the strong-coupling limit A — > 0, we find instead 

(-S.C.) = {- ^ (41) 
where Jo again is the closest multiple of 4 to J . 



2.6 Effective field equation 

An effective field equation can be obtained by substituting eq. (^91) into 
eq. (]40|), leading to 

"^p = ^E"^p' - -"sinM ^mp, + A/p J (42) 
pp' Vpp' / 



where we dropped the expectation value brackets for m^. This is equivalent 
to 



^L'TTip + Np = —Jp = —4a sin 



TT 



- (A^mp + Arup - Np) 



(43) 



We recognize on the left-hand side Stokes' law (plQl). Consequently, the new 
term on the right-hand side represents Jp, the density of monopole currents. 
It depends nonlinearly on itself and on the plaquette variable nip. In actual 
calculations, we use the numerical result from Fig. |I| in place of the sine. 
The periodicity of the right-hand side with respect to nip assures the residual 
gauge invariance under a global integer shift. 

To get a first impression of the meaning of eq. (|45|), we linearize the sine. 
This gives after some algebra 



nij, + AL = . 



1 - 27ra. 

Using Stokes' law (|l^), this is in continuum notation 



-J 



1 - 27ra 



-m = 



By taking the curl, one obtains 

1 - 27ra 



E 



Svra 



(44) 



(45) 



(46) 



which is a London relation between the electric field and the monopole current 
(cf. eq. (4) and (5) of |]l2l and eq. (15) of fl^; note also that the Dirac string 
on the right-hand side comes in naturally in our formulation). The London 
penetration length is in this approximation 



'1 - 2na 



(47) 



It duly diverges as a — >■ (weak-coupling limit). At a critical value of 



Ocrit = — 
ZTT 



(48) 
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the London penetration length becomes zero and then imaginary, thus sig- 
nalhng the breakdown of the hnearizing approximation. At this point, 
eq. P5| ) reduces to 

E = N , (49) 

i.e. the electric field is the Dirac string. 

In the strong-coupling limit, we can use without the linearization eq. (P|) 
for the right-hand side. This leads to 



- A^mp + A/p = -(A^mp + 4mp - A/p) + 4 
or equivalently 



Az.-Ap + mp 1 
4 2 



A ^mp — Ap + 4mp 1 

4 ^ 2, 



(50) 



(51) 



The brackets again indicate rounding down to an integer. This equation is 
in particular satisfied for mp = 0, which again corresponds to 



E = N 



(52) 



i.e. the Dirac string. 



3 Numerical solution and Discussion 

The nonlinear partial differential equation (|43|) can be solved numerically 
in certain circumstances. We used a simple iterative method in which the 
equation is solved by simultaneous over-relaxation for a fixed right-hand side. 
The right-hand side is then adjusted to the new values of {nip), and the 
process repeated until convergence. We found that this method converges 
provided that the lateral extension of the lattice (i.e. the y dimension) is 
not too large; otherwise spontaneous oscillations set in. In particular, the 
difficult part here is the strong-coupling region that is highly nonlinear. 

Fig. ^|-^ show the electric field calculated from eq. (^) for different val- 
ues of A at a constant string length of 40 lattice units. The lattice size was 
between 60x100 at weak coupling and 60x22 at strongest coupling. The solu- 
tion varies between the string-like strong-coupling limit and the Coulomb-like 
weak-coupling limit. 
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Figure 2: Electric field across the string axis for different values of A 
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Figure 3: Electric field along the string axis for different values of A 
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The actual density of monopole charges in a cut through the string is 
shown in fig. ^ Their concentration is largest just next to the string axis 
where the field is confined, and then drops sharply. This justifies the inter- 
pretation that the monopole charges create a "wall" that confines the electric 
field. 

In fig. H, we have calculated the root-mean-square radius 



(where y is the distance from the string axis) orthogonal to the string axis of 
the fiux tube as a function of A and of the position along the string axis. For 
large A, it is constant along most of the string axis, indicating the formation 
of a fiux tube. For small A, it smoothly reaches a maximum halfway between 
the charges as a result of the 1/r-law. 

Fig. ^ shows the root-mean-square radius for different lengths of the 
string. We do not see the Inr dependence predicted from studying 



quantum fiuctuations of the string. The width of the string stays constant, 
as would be expected from strong-coupling expansions. This is an indication 
that quantum effects will modify eq. (|43|) . 

The numerical solution of eq. (^) exhibits flux tube formation as is ex- 
pected for conflnement. As A is increased and the weak-coupling limit is 
approached, the flux tube disappears. On the other hand, it is expected 
from studies of the effective monopole action that 2 + 1-dimensional U{1) 
lattice gauge theory should not exhibit a phase transition to deconflnement 



2^ . Both results can be reconciled when considering that conflnement of the 
electric fleld does not at flrst hand prove conflnement of charges or provide 
for the calculation of the string tension. For the latter, it has to be shown 
that the ground-state energy of the Hamiltonian grows proportional to the 
distance of the charges. However, when truncating to single-plaquette states, 
any information about the correlation energy is lost, and the ground-state 
energy cannot be calculated rehably (cf. eq. (P?])). For such a calculation, 
blocked ground states as used in real-space renormalization must be utilized, 
and the renormalization group, which is supposed to be the prime cause for 
conflnement of charges, is brought back in. Quantum effects manifested by 
a roughening transition might then modify eq. (^) to reinforce conflnement 
of the electric fleld. 
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Figure 6: Electrical field strength \E\ at A = 8 




Figure 7: Monopole density in a cut across the string axis for different A 
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Figure 8: Root-mcan-square radius of the string along the string axis for 
different values of A 




Figure 9: Root- mean-square radius of the string along the string axis for 
different distances of the charges 
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4 Conclusions 



We have shown that a semiclassical confining field equation can be directly 
derived from the Hamiltonian formulation of lattice gauge theory. Its lin- 
earized form is a London relation between the electric field and the monopole 
current. Numerical calculations show the formation of a fiux tube for the elec- 
tric field between two charges. The wall of the fiux tube is formed by effective 
monopole currents that result from the nonlinearities of the field equation. 

These results have been obtained using a rather crude discorrelated vac- 
uum state. One of the major questions that accompanies this approach is 
how this equation is modified by correlations, i.e. quantum effects, and how 



this is related to a possible roughening transition ||26[. The validity of the 
London relation for strong coupling has been demonstrated in Monte Carlo 
calculations. For weak coupling, renormalization group analysis predicts con- 
tinued confinement of charges while the semiclassical field equations converge 
to the limit of classical electrodynamics. Here the inclusion of correlations 
(and consequently the renormalization group behavior) is needed to calcu- 
late the string tension and show confinement of charges. The most promising 
way to do this is to extend the truncation scheme to correlated blocks and 
perform a Hamiltonian real-space renormalization group analysis. 

A direct check of the validity of eq. (^31) is possible by comparing the 
results to Monte Carlo simulations. In particular, we are investigating the 



guided random-walk ground-state ensemble projector method [|2^ that is 
based on the Hamiltonian approach. 

Further extension of the model is to 3 -|- 1 dimensions and to the SU (2) 
gauge group. In the extension to 3 + 1 dimensions, one encounters constraints 
in the plaquette representation (Bianchi identities). While this will probably 
affect correlations, it remains to be seen whether the effective field equations 
changes qualitatively, and how this can be related to the maximally Abelian 
gauge which has been used in some of the Monte Carlo calculations. 

Finally the extension to finite temperature can be considered to study 
the deconfinement transition. 

This work was supported by Deutsche Forschungsgemeinschaft (DFG). 
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